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Abstract 

The role of chiral transformations in effective theories modeling Quantum Chromo Dy- 
namics is reviewed. In the context of the Nambu-Jona-Lasinio model the hidden gauge 
and massive Yang-Mills approaches to vector mesons are demonstrated to be linked by a 
special chiral transformation which removes the chiral field from the scalar-pseudoscalar 
sector. The chirally rotated axial vector meson field (A^) transforms homogeneously un- 
der flavor rotations and may thus be dropped without violating chiral symmetry. The 
fermion determinant for static meson field configurations is computed by summing the 
discretized eigenvalues of the Dirac Hamiltonian. It is discussed how the local chiral 
transformation loses its unitary character in a finite model space. This technical issue 
proves to be crucial for the construction of the soliton within the Nambu-Jona-Lasinio 
model when the chirally rotated axial vector field is neglected. In the background of this 
soliton the valence quark is strongly bound, and its eigenenergy turns out to be negative. 
This important physical property which is usually generated only by non-vanishing axial 
vector is thus carried over by the simplification = 0. 
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1. Introduction 



Since a solution to Quantum Chromo Dynamics (QCD) is not yet available one has 
to recede on models in order to explore processes described by the strong interaction. 
These models are usually constructed under the requirement that the symmetries of 
the underlying theory, i.e. QCD are maintained. In this context chiral symmetry and 
its spontaneous breaking are of special importance. In this article we will explore a 
special chiral transformation when topologically non-trivial meson field configurations 
like solitons are involved. To begin with, let us briefly review the relevance of chiral 
symmetry on the one side and solitonic field configurations on the other in the context 
of strong interactions. 

QCD can be extended from SU{3) to SU{Nc) where Nc denotes the number of 
color degrees of freedom. It was observed by 't Hooft [|I] that in the limit Nc oo 
QCD is equivalent to an effective theory of weakly interacting mesons. Subsequently 
Witten 0] conjectured that baryons emerge as solitons of the meson fields within this 
effective theory. Stimulated by Witten's conjecture much interest has been devoted to 
the description of baryons as chiral solitons during the past decade |P, ^ . In the soliton 
description of baryons the chiral field and in particular its topological structure play a 
key role. The topological character of the chiral field especially endows the soliton with 
baryonic properties like baryon charge and spin [§. This comes about via the chiral 
anomaly which is a unique feature of all quantum field theories where fermions live in 
a gauge group and couple to a chiral field. The fermionic part of such a theory has the 
generic structure 



Zp[^] = J D-^D-^ exp 



= Det {ip - mo) (1.1) 

where rho denotes the current quark mass matrix which will be ignored in the ongoing 
discussions. Furthermore 

ip=i^-^ = i{^ + D- MPr- m'^Pl (1.2) 

represents the Dirac-operator of the fermions in the external Bose-field $. $ in general 
contains vector V^, axial vector field^ as well as scalar S and pseudo-scalar fields P 

r, = V^ + A^j,, M = S + zP = d^U- (1-3) 

Here Pr/l = |(1 ± 75) are the chiral projectors. Accordingly one defines left(L)- and 
right (R)-handed quark fields: '^l,r = Pl,r^ ■ The chiral field U is defined via the polar 
decomposition of the meson fields 

t/ = &ij = exp(ze). (1.4) 

The chiral anomaly arises because there is no regularization scheme which simulta- 
neously preserves local vector and axial vector (chiral) symmetries. In renormalizable 



°In our notation and are anti-hermitian. 
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theories the chiral anomaly can be calculated in closed form and is given by the Wess- 
Zumino action 0, |^ 

^=Swz = ^J {Udu'ff (1.5) 

where the notation of alternating differential forms has been used. denotes a five 
dimensional manifold whose boundary is Minkowski space. Obviously the chiral anomaly 
is tightly related to the chiral field since the Wess-Zumino action vanishes when the chiral 
field disappears {U = 1). The chiral anomaly is, however, not merely a technical artifact 
but has well established physical consequences. In the meson sector it gives rise to the 
so-called "anomalous decay processes" like e.g. vr — > 27 and uj — > Sir. In the soliton 
sector the chiral anomaly requires for A'",^ = 3 the soliton to be quantized as a fermion 
and endows the soliton with half integer spin and integer baryon numberti J5|, ^. 

For many purposes it is convenient to perform a chiral rotation of the fermions |10, 11 

^ = n^ with n = PL^L + PR^R, i.e. ^l,b = ^l,r'^l,r. (1-6) 

This transformation defines a chirally rotated Dirac-operator 

^ip\i/ = ^ip^ (1.7) 

which acquires the form 

ip= nUpn^ = i-f^ (9^^ + W + - s. (1.8) 

The chiral rotation has removed the chiral field from the scalar pseudo-scalar sector of 
the rotated Dirac operator ip. As a consequence the vector and axial vector fields now 
become chirally rotated 



V,-A, = ^L{d, + V,-A,)d. (1.9) 

Even in the absence of vector and axial vector fields in the original Dirac operator (V^ = 
Afj, = 0) the chiral rotation induces vector and axial vector fields 

MV, = A, = 0) = v, = 1 (^eH^^ei + a^^d 

A,iV, = A, = 0) = a, = U^j,d,^l-^^d,d)- (1-10) 



For the soliton description of baryons the chiral field is usually assumed to be of the 
hedgehog type 

U = exp {ie{r)T ■ r) (1.11) 

The non-trivial topological structure of this configuration is then exhibited by the bound- 
ary conditions 0(0) = — nvr and O(cxd) = 0. The chiral field thus represents a mapping 



''For this proof it is mandatory to consider flavor SU(3). 
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from the compactified coordinate space (all points at spatial infinity are identified) to 
SU(2) flavor space, i.e. S^. The associated homotopy group, 113(5'^), is isomorphic 

to Z, the group of integer numbers. The isomorphism is given by the winding number 

(0(0) — 0(oo)) /tt = —n. Assuming the unitary gauge (^| = ^r) the induced vector field 



is of the Wu-Yang form [11 



t;o = 0, v, = ivf-^, y- = e'kaf^^ (1.12) 

with the profile function G{r) given by the chiral angle G(r) 

G{r) = -2sin2^. (1.13) 

For odd n the topological non-trivial character of the chiral rotation is also reflected by 
a non- vanishing value of the induced vector field at r = 0. The induced axial vector field 
tti = ialr"" /2 becomes 

n - - (r^t/ \ sin0(r)\ ^ sin9(r) ,^ 
^ r y r 

where the prime indicates the derivative with respect to the argument. 

The use of the chirally rotated fermions is advantageous since the rotated Dirac- 
operator does no longer contain a chiral field and its determinant is hence anomaly free. 
The chiral anomaly, however, has not been lost by the chiral rotation but is hidden in 
the integration measure over the fermion fields. In fact, as observed by Fujikawa |T^, a 



chiral rotation of the fermion fields gives rise to a non-trivial Jacobian of the integration 
measure 

D^D^! = J{U)D^D^ (1.15) 
which is precisely given by the chiral anomaly 

.J{U) = exp (iA) (1.16) 

Thus we have the relation 

-iTi log ip = -iTi log ip+ A (1.17) 

In many cases it is convenient to work with the chirally rotated fermion fields because of 
the absence of the anomaly from the fermion determinant. 

As already mentioned above the chiral anomaly or equivalently the non-trivial Jaco- 
bian in the fermionic integration measure arise due to the need for regularization, which 
introduces a finite cut-off. In the regularized theory the chiral anomaly can be evaluated 
in a gradient expansion. In leading order the anomaly is then given by the Wess-Zumino 
action ( [1.51 ). There are, however, higher order terms which are suppressed by inverse 
powers of the cut-off. In renormalizable theories where the cut-off goes to infinity, these 
higher order terms disappear and the chiral anomaly is known in closed form. In non- 
renormalizable effective theories, however, the cut-off of the regularization scheme has 
to be kept finite and acquires a physical meaning, indicating the range of validity of the 
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effective theory. In this case the higher order terms of the gradient expansion do not 
disappear but contribute to the anomaly which is then no longer available in closed form. 

Furthermore, when the soliton sector of such effective mesonic theories is studied it is 
not sufficient to only consider the leading and sub-leading contributions from the gradient 
expansion and one has to perform a full non-perturbative evaluation of the fermion 



determinant [13|. The non-perturbative calculations have to be performed numerically 



14], with the continuous space being discretized. This is the case for both, coordinate 



and momentum space. Also in the non-perturbative studies of the soliton sector of the 



effective theory the use of the chiral rotation is in many cases advantageous |Tl|. As 
noticed above for the soliton description of baryons the topological nature of the chiral 
field is crucial. Actually, topology is a property of continuous spaces (manifolds) and it 
is a priori not clear whether the chiral rotation with a topologically non-trivial chiral 
field can be represented in a finite dimensional and discretized model space used in the 
numerical calculations. In this respect let us recall that a single point defect in a manifold 
changes its topological properties already drastically. 

The purpose of this paper is twofold. First we wish to present a study of the chiral 
rotation in non-perturbative soliton calculations where the fermion determinant has to 
be numerically evaluated in the background of a topologically non-trivial chiral field. It 
will be demonstrated how this rotation influences the choice of the boundary conditions 
for the eigenstates of the Dirac Hamiltonian. Second we will make use of the fact that 
in the chirally rotated formulation the axial vector field may be set to zero {i.e. = 0) 
without spoiling chiral symmetry. The resulting soliton configuration is constructed and 
compared with various solutions in the unrotated formulation. For definiteness we will 
use the Nambu-Jona-Lasinio (NJL) ||15| model as a microscopic fermion theory which 
shares all the relevant properties of chiral dynamics with QCD. Its bosonized version 
gives a quite satisfactory description of mesons [|10] and also of baryons when the soliton 
picture is assumed 



The organization of the paper is as follows: After these introductory remarks we 
review the importance of the local chiral rotation ( |1.6| ) for the extraction of meson prop- 
erties from the NJL action. In that section also the chirally invariant elimination of the 
axial vector meson is described. In section 3 we review the computation of the soliton 
solution to the NJL model of pseudoscalar mesons. This may straightforwardly be gen- 
eralized to the inclusion of other mesons as long as the Euclidean Dirac Hamiltonian 
remains Hermitian. Section 4 is devoted to the study of the local chiral rotation for topo- 
logically non-trivial field configurations and its influence on the soliton solution. The 
computation of the soliton solution with the rotated axial vector meson field being elim- 
inated is described in section 5. Previously it has been shown that the presence of the 
axial vector field is responsible for a strong binding of the valence quark and the valence 



quarks' eigenenergy being negative ^7\, [T^. We will in particular examine whether this 
important piece of information is maintained after eliminating A^. This is not obvious 
because simply setting A^ = leads to a positive valence quark energy for a reasonable 
choice of parameters A concluding discussion is given in section 6. Some technical 
remarks on the boundary conditions for the Dirac spinors are left as an appendix. 

2. Properties of local chiral transformations 

As indicated in the introduction the chirally rotated formulation of the NJL model 
is suited to investigate properties of (axial-) vector mesons. In the present section we 
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will therefore briefly review the results and illuminate the connection with the hidden 
gauge symmetry (HGS) and massive Yang-Mills (MYM) approaches for the description 
of (axial-) vector mesons. Many of the results reported in this section are taken from 
earlier works [l^, |ll], Nevertheless we repeat these results here in order to put 
our work into perspective and have the paper self-contained. 

In the chirally rotated formulation the bosonized version of the NJL model action 
reads 



A 



NJL 



-zTr \ogip+ A — 



1 



4Gi 



tr - mn - 



1 



AG, 



-tr 



[2.1] 



with the constituent quark mass m being the vacuum expectation value of the scalar 
field S i.e. (S) = m. Again we have discarded terms proportional to the current quark 
mass matrix. The chirally transformed (axial-) vector fields are defined in eqns ( |1.9D a nd 
( |1.10| ). Next we have to face the fact that the functional trace of the logarithm in ( |2.1|) is 
ultra-violet divergent and thus needs regularization. This is achieved by first continuing 
to Euclidean space (xq = —ixi) and then representing the real part of the Euclidean 
action by a parameter integral 



^Trlog to^' 



ds 



2 Jl/A2 



-exp 



(2.2) 



which introduces the cut-off A. This substitution is an identity up to an irrelevant 
additive constant for A oo. The Euclidean Dirac operator is obtained by analytical 



continuation of pto Euclidean space. The prescription ( |2.2| ) is known as the proper-time 
regularization For the purpose of the present paper it is sufficient to only consider 
the normal parts of the action. We may therefore neglect the imaginary part as well as 
the anomaly A. Thus the actual starting point of our considerations is represented by 



^NJL = -^/J^.TTrexpf-4l^ 



1 



4G9 



-tr 



. 2 



tr S - 



Afi — dfj, 



(2.3) 



In order to extract information about the properties of the (axial-) vector mesons com- 
monly a (covariant) derivative expansion of the fermion determinant is performed. We 
choose to consider a covariant derivative expansion since, in contrast to an on-shell de- 
termination of the parameters [E^], it preserves gauge invariance and does not lead to 



artificial mass terms. Furthermore, this procedure leaves the extraction of the axial- 
vector meson mass unique. Continuing back to Minkowski space and substituting the 



scalar field S by its vacuum expectation value yields the leading terms EO, 111 



NJL 



2gl 



tr 



+ A^ 



tr 



+ ... .(2.4) 



Here V^y and A^^ denote the field strength tensors 



V, 



A 



d^A,-d,A, + [%,A,] + [A^,V,] 



(2.5) 
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of chirally rotated vector and axial-vector fields, respectively. Obviously in our convention 
these fields contain the coupling constant gy which in the proper-time regularization is 
given by 



9v 



2N, 




(2.6) 



For the description of the pion fields we adopt the unitary gauge for the chiral fields: 

= ^R- The pions come into the game by the non-linear realization ^ = exp (zr ■ tt//). 
Then the last term in eqn (|2.4| ) contains the axial-vector pion mixing which is eliminated 



by a corresponding shift in the axial field: A'^ = + {igyW? /12fG2)d^T -tt. This 

shift obviously provides an additional kinetic term for the pions and thus effects the pion 
decay constant 



f 

J TT 



with the (axial-) vector meson masses 



Ml 



Ml 



4:MiG2 



Ml 



AG. 



and 



Ml 



Ml + 6m2. 



(2.7) 



(2.8) 



This brief summary of known results has demonstrated the usefulness of the chiral ro- 
tation ( [I.6| ) especially in the context of the derivative expansion since it eliminates the 
derivative of the chiral field from the fermion determinant ( |1.8|) . 

The Lagrangian of the hidden gauge approach can be obtained from (|2.4|) by the 
following approximation. One neglects the kinetic parts for the axial-vector field A'^ 
which leaves this field only as an auxiliary field. This allows to employ the corresponding 
equation of motion to eliminate A'^ resulting in 



1 



2^7 



2 tr v;t 

V 



a flu {y,-v,)^ --fliia 



In the work of Bando et al. |^ a was left as an undetermined parameter, 
fixed in terms of physical quantities 

Ml 



Ml 



(2.9) 



Here it is 



(2.10) 



Assuming the constituent quark mass m = My / \^ not only yields the Weinberg relation 
Ma = V2Mv [H but also the KSRF relation a = 2 

Alternatively one might apply the same manipulations to the formulation in terms of 
the unrotated fields ( |1.3|) . Then the chiral field still appears in the fermion determinant 
and one has to deal with the covariant derivative 



V^U = d^U+[V^,U]-{A^,U} 



(2.11) 



The leading terms in the Lagrangian can readily be obtained [T^ 



C 



3m2 
2^ 



tr (vMV^'u'f 



tr V. 



A 



AG. 



-tr 



(2.12) 
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which exactly represent the massive Yang-Mills Lagrangian p, |27] 
(axial-) vector fields according to ( pTOl) and noting that |^| 



tr Al = -^ti ( V„,UV^U^ 



Transforming the 



(2.13) 



one immediately observes that (|2.12|) and ( p.4|) describe the same physics. In the language 
of the NJL model the identity of the HGS and MYM approaches stems from the invariance 
of the module of the fermion determinant under the special chiral rotation 



1.6) 



After we have seen how the equivalence of the hidden gauge and massive Yang Mills 
Lagrangians emerge from the NJL model we next explore the transformation properties 
of the fields under fiavor rotations gi, Qr- These properties ploy a key rolejn order to 
demonstrate that the elimination of the axial vector meson field by setting = does 
not violate chiral symmetry. The fiavor rotations are defined for the unrotated left- and 
right-handed quark fields 



^1 



and 



R 



Qr'^r- 



(2.14) 



The term which describes the coupling of the quarks to the scalar and pseudoscalar 
mesons is left invariant by demanding 



9Ld^U9R 



(2.15) 



which introduces the hidden gauge transformation h |10 



l9l 



and 



[2.16) 



Obviously the scalar fields transform homogeneously under the hidden gauge transfor- 
mation. In this context it is important to note that h may not be chosen independently 
but rather depends on the gauge adopted for the chiral fields. Consider e.g. the unitary 

gauge d = C,R = C,- This requires the transformation property 



For vector type transformations gi = gR 



(2.17) 



gv this equation is obviously solved by h 
Contrary, for axial type transformations gi = g\i = gA h is obtained as the solution to 
gA^h = h^^gA which depends on the field configuration ^. Thus even for global fiavor 
transformations gA,v the hidden gauge transformation h may be coordinate-dependent 
for coordinate dependent ^(x) [| 



29[| . The unrotated (axial-) vector fields are required 



to transform inhomogeneously under the fiavor rotations 



gR {d^ + V^ + A^) and 



gL{d^ + V^-A^)gl. (2.18 



It is then straightforward to verify that the fiavor transformation of the rotated fields 



only involves the hidden symmetry transformation h 



L,R 



L,R 



and 



^ /iti^/i. (2.19) 
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The fact that transforms homogeneously has the important consequence that (as 
already mentioned) one can set = without breaking chiral symmetry b cause the last 
relation in ( p.l9| ) does not induce any inhomogeneity. Furthermore the vector mesons 
are not affected by this choice. Let us next examine the NJL model defined by = 
in more detail. The corresponding Dirac operator reads 



and the derivative expansion (p.4|) simplifies to 



^2.20) 



-tr V.i - 



AG- 



-tr 



+ a. 



+ 



(2.21) 



Identifying with the physical p-meson field determines the coupling constant G2 

27r2 



9l 



4G2 G^V (0, f ) 



(2.22) 



for = 3. One may rewrite the last term in equation ( |2.2lD in terms of the chiral field 
U (see also eqn. (|2.13|) ) 



tr a, 



--tr d^Ud^U 



(2.23) 



which provides an additional relation for G2 in terms of the pion decay constant 

1 



f 

J IT 



This finally implies 



4G2 

.2 / A 2\ 



(2.24) 

It should be remarked that the above 



^ 87rV.Vr(0,mVA 
derived relations are also obtained when (in a derivative expansion) is set to zero and 
Vfj_ identified with the p-meson. This, however, is not surprising since these relations for 
the vector mesons stem from the part of the Lagrangian which does not contain pion 
fields. In the absence the pion fields and are identical, cf. eqn ([1.91). 
The derivative expansion (p.21|) also contains the pnir vertex 



QpTTTT 



-tr 



^2"^^ ' ' 2G2 
Expanding the RHS in terms of the physical fields and tt yields 



rrir 



gpnn 



V2U 



[2.25) 



[2.26) 



which falls short a factor ^/2 of the KSRF relation. This seems to indicate that some 
relevant information is lost by setting = 0. Whether this is also the case in the soliton 
sector will be explored in section 5. 

In phenomenological vector meson models frequently a term like 



tr 



[2.27) 



9 



is added p9|, to reproduce the empirical value for (7^7^77 when setting = 0. This 
term is invariant under global chiral rotations. The absence of this term in the NJL 
model causes the improper prediction for this coupling constant. 

3. The NJL soliton 

In the baryon number one sector the NJL model has the celebrated feature to possess 
localized static solutions with finite energy, i.e. solitons |13, 14|. Here we wish to briefly 
review this solution for the case of pseudoscalar fields. 

For static field configurations it is convenient to introduce a Dirac Hamiltonian 

n = a-p + P (pRi{m + PLiHT.)A (3.1) 



where we have assumed the unitary gauge (i.e. ^1 = ^i? = 0- This Hamiltonian enters 
the Euclidean Dirac operator via 

i^pj^ = - n. (3.2) 

For static mesonic background fields the fermion determinant can conveniently be ex- 
pressed in terms of the eigenvalues of the Dirac HamiltonianB 

n-^^ = e^^^. (3.3) 

These eigenvalues are functionals of the mesonic background fields. Depending of the 
specific boundary condition (which fixes the quantum reference state) to the fermion 
fields in the functional integral (|1 . 1| ) , the fermion determinant (|1.1|) contains in general 
besides a vacuum part also a valence quark part A^^^ [|13| 

^ = ^0 + (3.4) 

The valence quark part arising from the explicit occupation of the valence quark levels 
is given by 

^val ^ _^val[^]y ^ ^val[^] = NcJ2 V,K\. (3.5) 

Here r/^ = 0, 1 denote the occupation numbers of the valence (anti-) quark states. These 
have to be adjusted such the total baryon number 

B = T.{%- ^sgn (e^)) (3.6) 

equals unity. The vacuum part is conveniently evaluated for infinite Euclidean times 
(T — >■ oo) which fixes the vacuum state as the quantum reference state (no valence quark 
orbit occupied). For the present considerations it will be sufficient to evaluate the real 
vacuum part 

A'n = ^Tr log^l^i^- (3.7) 



"The treatment generalizes to more complicated field configurations as long as Ti. is Hermitian (c/. 
section 5). 
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Since for static configurations one has [dr,'H]=0 and thus IJ^ePe = —d"^ + 'H^- Then it 
is straightforward to evaluate the real part of the fermion determinant in proper-time 
regularizationi 



1/A2 

Nc r°° dz ^ r°° ds 



-T^ -exp + el]). (3i 



s 

The temporal part of the trace has become the z integration. As this integral is Gaussian 
it can readily be carried out yielding 

This expression allows to read off the static energy functional E[C,] since ~^ ^TE^[^] 
as T — > cx) 

E'[^] = ^r ^Eexp(-4). (3.10) 



Finally the total energy functional is given by 

E[i] = E^^\i\ + E'[i\ - E'[i = I] (3.11) 

which is normalized to the energy of the vacuum configuration ^ = 1. In the chiral limit 
(m^r = 0), which we have adopted here, the meson part of the action does not contribute 
to the soliton energy. The chiral soliton is the ^ configuration which minimizes E[^] and 
the minimal E[^] is then identified as the soliton mass. 

To be specific we employ the hedgehog ansatz for the chiral field 

e(r) = exp Qr • f e(r)) (3.12) 

while the scalar fields are constrained to the chiral circle, i.e. (S) = m. Substituting this 
ansatz into the Dirac Hamiltonian (|3.1|) yields 



7i = CK • p + /3m (cos9(r) + i75T ■ f sin6(r)) . (3.13) 

The stationary condition 6E[^]/6^ = is made explicit by functionally differentiating 
the energy-eigenvalues with respect to 9 



— m 



6Q{r) 



J dVL ^J(r)/3 (-sine(r) + i-f^r ■ r cose(r)) ^^(r). (3.14) 



This leads to the equation of motion [|T^ 

cos6(r) tr J dil ps{r,r)i'y5T ■ r = sin9(r) tr J dil ps{r,r) (3.15) 

''The imaginary part does not contribute for the field configurations under consideration because H 
is assumed to be Hermitian. 
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where the traces are over flavor and Dirac indices only. According to the sum (|3.11|) the 
scalar quark density matrix ps{x, y)) = {q{x)q{y)) is decomposed into valence quark and 
Dirac sea parts: 

psix,y) = pj\x,y)+pY{x,y) 
pfix,y) = Y.'^^{x)rj^^^{y)sgn{e^) 



PTi 



x^y) = ^E^MNerfc(|^|)^^(?/)sgn(e^). (3.16) 



The explicit form of the eigenfunctions "^fxix) as well as remarks on the appropriate 
boundary conditions may be found in the appendix. 

4- The chirally rotated fermion determinant 

In the appendix it is demonstrated that the normalizable solutions to the free Dirac 
equation with spherical boundary conditions represent a pertinent basis for the diagonal- 
ization of the Dirac Hamiltonian with the soliton present. This property is based on the 
fact that the Hamiltonian (|3.13|) is free of singularities. In this section we will explain 
how singularities appearing in a Dirac Hamiltonian influence the choice of basis states. 
Let us for this purpose consider the Hamiltonian for the chirally rotated quark fields 
^ = Q{e)^ {cf. eqns (|LT2D-(|LT4D and ref. 



nR = ni<d)nnUQ) = (X-p + Pm--((T-r)(T-r)(e'(r)--sme(r 

2 V r ' 

- — (a- ■ T)sine(r) --a-(r x r) sin^ [ ) (4.1) 
2r r \ 2 J 

since in unitary gauge f2(9) = cos(6/2) +ij^T-r sin(6/2). Obviously the 6-dependence 
in the Hamiltonian has been transferred to induced (axial-) vector meson fields. As ex- 
pected the rotated Hamiltonian, Tin, contains an explicit singularity in the 1/ra ■ (f x 
t) sin^ (^2^) term at r = 0. Additionally there are "coordinate singularities" in the 
expressions involving r. All these singularities appear because the "coordinate singular- 
ity" in f2(G) is not compensated by corresponding values of the chiral angle B. Stated 
otherwise: the transformation Q{Q) with 0(0) — 0(cxo) = —nir is topologically distinct 
from the unit transformation. Although Q{<d) represents a unitary transformation it is 
then not astonishing that a numerical diagonalization 

Hr^^ = e^^^ (4.2) 

in the basis of the free Hamiltonian does not render the eigenvalues of the unrotated 
Hamiltonian, Ti, {i.e. 7^ e^) despite the relevant matrix elements being finite. This 
finiteness is merely due to the factor in the volume element. One might suspect 
that the Hamiltonian Ti.2R = obtained by a 2G rotation has the same 

spectrum as Ti, since 7^2/? is free of singularities. Although this behavior is exhibited by 
the numerical solution for the low-lying energy eigenvalues, the topological character of 
the transformation has drastic consequences for the states at the lower and upper ends 
of the spectrum in momentum space. Adopting the same basis states for diagonalizing 
H and 'H2R one observes that the eigenvalues of Ti.2R are shifted against those of H, 
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i.e. the most negative energy eigenvalue is missing while an additional one has popped 
up at the upper end of the spectrum. Up to numerical uncertainties the eigenvalues in 
the intermediate region agree for both Ti. and Ti.2R. This behavior is sketched in figure 
4.1 and repeats itself for H4R = Thus the chiral rotation represents 



another example of the so-called "infinite hotel story" |3ll] which is an interesting feature 
reflecting the topological character of this transformation. 

Let us now return to the problem of diagonalizing Tin and restrict ourselves for the 
moment to the channel = 0^. The grand spin operator G is defined in the appendix 
( |A.lj ). It should be stressed that [f^, G] = 0. Thus the local chiral rotation may be 
investigated in each grand spin channel separately. 

At r = the chiral rotation 

n(r = 0) = -i (r • f) 75 (4.3) 

obviously exchanges upper and lower components of Dirac spinorsi. The corresponding 
wave-functions \^|f'"'")(r = 0) = ^l{r = 0)\E'|f'"^-'(r = 0) (see eqn ([A.9|)) can obviously not 



be represented by the free basis as in general the lower component of tilde'^^^'~^\r = 0) 
is different from zero. However, the lower components of the eigenstates of the free 
Hamiltonian in the = 0^ channel have this property [cf. eqns ( |A.3| , |A.5D ). It is thus 



clear that the local chiral rotation is not unitary in a finite model space. Furthermore 
the equality = cannot be gained without additional manipulations. Of course, this 
"non-unit arity" is completely due to the topological character of this rotation. This 
problem can be avoided by defining a basis in the topologically non-trivial sector via 

^,0 = ^^(0)^mo (4.4) 

with being the solutions to the free unrotated Hamiltonian. (p represents an auxiliary 
radial function satisfying the boundary conditions 0(0) = —it and 0(-D) = 0. E.g. we 
may takeS 



0(r) = — TT 1^1 — — j exp {—tmr) (4.5) 

with t being a free parameter. The diagonalization of Tin in the basis \E'o is equivalent 
to diagonalizing 



1 

+2 



a ■ p + m(3 (cos 0(r) + i75T ■ r sin 0(r)) 

cr ■ TT ■ r 



0'(r) - e'(r) + - sin(e(r) - 0(r)) 



+— sin(0(r) — Qir))cr ■ r — [1 — cos(0(r) — Mr))] a ■ (f x r) (4.6) 

2r 2r 

in the standard basis {\E'^o}- At this point it should be stressed again that cf) is not 
a dynamical field but merely an auxiliary field which transforms the basis such as to 
eliminate the 1/r-singularities from the Dirac Hamiltonian. This property is completely 



°In the standard representation 75 



1 

1 



denotes the size of the spherical cavity which serves to discretize the momentum eigenstates, cf. 
the appendix. 
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determined by the boundary values of 0. Thus the results ought to be independent of 
the parameter t. We have confirmed this property numerically. The total energy of the 
soliton varies only by fractions of MeV in a wide range of t. This is, of course, negligibly 
small since the inherited mass scale is given by the constituent quark mass, m which is of 
the order of several hundred MeV. By using the locally transformed basis (|4.4| ) also the 
wave-functions corresponding to the eigenstates of Hr agree reasonably well with the 
rotated wave-functions ^2(0)^^^ of the original Dirac Hamiltonian Ti.. It should be noted 
that a large number of momentum states is required to numerically gain this result. This 
is not surprising since in order to represent the functional unity an infinite number of 
momentum states is needed. 

At this point we want to add a remark on a further application of the chiral rotation. 
Using = ±7r, i.e. a chiral transformation Q = ±z75T ■ f, the rotated Hamiltonian ( |4.1| ) 
is given by 

n^^^"" = ct ■ p + pm - ■ {r X t) (4.7) 



whereas the unrotated one (see eqn (|3.13|) ) is simply 7-^®=='='^ = a ■ p — jSm. The latter is 
diagonalized straightforwardly. It has the same eigenvalues as the free Hamiltonian and 
the eigenf unctions differ by the substitution m — >■ —m. Therefore the Hamiltonian 

( [4.7D has also the eigenvalues of the free Hamiltonian. The eigenf unctions are given by 

^®=±'^ = l]t(e = ±71)-^^=^^ = ±ij5T ■ f^®=±^. I.e. the eigenf unctions of (|3) are 
easily constructed: start with the free eigenfunctions ( |A.3 - A.5 ), substitute m —>■ —m and 



apply ±i75T ■ r. This constitutes just one example how the chiral rotation ( |1.6| ) can be 
used to diagonalize non-trivial operators. 

Before turning to the detailed discussion of the equation of motion in the rotated 
system we would like to mention that the self-consistent profile being obtained from the 
unrotated problem also minimizes the soliton mass in the chirally rotated frame. Stated 
otherwise: each change in this profile function leads to an increase of the energy obtained 
from the eigenvalues of Hr. 

As in the formulation within the unrotated frame the equation of motion is gained by 
extremizing the energy functional ( 3.11 ). The energy eigenvalues in the rotated frame 



however, exhibit a different functional dependence on the chiral field. The functional 
derivative of these eigenvalues with respect to the chiral angle reads 



5Q{r) 



/I C7 1 ~ + ~ 

dn \-(— + -cos&(r))r'^m(r)(T ■ rr ■ f^„(r) (4.^ 

-cos6(r)^|(r)cr ■ r^^(r) - -sin6(r)^f J(r)Q; ■ (f x r)^^(r)| 



with ^ ^[t^ being the eigenstates of the rotated Dirac Hamiltonian ( ^4.1| ). The derivatives 
( [4.8[ ) enter the stationary condition for the energy functional resulting in the equation of 
motion 

Ai(r) + AT(r)cos0(r) - 1/(r)sin0(r) = 0. (4.9) 

In order to display the radial functions AljAx and V it is convenient to introduce the 
charge density pc = {Q{x)q{y)^) = p}f + Pc"^ involving the eigenstates of Hr (cf. eqn 
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A 



^!(2/)sgn(e^). 



(4.10) 



The radial functions Ai,At and V are of axial-(Ai r) and vector(V^) character 

Id r 

Aiir) = -— tr d^l r"^ pc(r,r)cr ■ TT ■ r (4-11) 
r or J 

Arir) = tr J dn pc{r, r) [cr ■ rr ■ r - cr ■ r] (4.12) 

V{r) = ti J dn pcir,r)a. ■ (f x r). (4.13) 

It should be noted that Al does not depend on the auxiliary field 0. It is then straight- 
forward to verify that for any meson configuration A^ and A^ satisfy the relation 



AL{r = 0) = -(-l)MT(r = 0) 



(4.14) 



where k is defined by the value of the auxiliary field (p at the origin 0(r = 0) = kn. The 
vector type radial function V vanishes at the origin. Thus the equation of motion ( ^.9| ) 
together with the relation ( [4.14| ) yield the boundary condition 0(r = 0) = (2n + l)7r for 
k = 1. This is stronger than the boundary condition derived from the original equation of 
motion ( |3.15D which also allows for even multiples of vr for G(r = 0) since tr / dQ ps{r = 
0, r = 0)275T ■ f = Assuming the Kahana-Ripka boundary conditions for the 
unrotated basis states \l/^o similar considerations for r = D show that Q{r = D) = 21tt 
since A^lr = D) = —A^ir = D) as long as (t){D) = 0. Obviously the topological charge 
associated with the chiral field in the hedgehog ansatz (0(r = 0) — 0(r = D))/tt can 
assume odd values only when k is odd. In particular = is prohibited in the case of unit 
baryon number. Thus the study of the boundary conditions in the chirally transformed 
system corroborates the conclusion drawn from investigating the eigenvalues and -states 
of Tin that it is mandatory to also transform the basis spinors and in particular the 
boundary conditions to the topological non-trivial sector. 

We would also like to remark that this kind of singularities does not only appear when 
the Dirac Hamiltonian is considered. Such topological defects have also caused problems 
in Skyrme type models when fluctuations off vector meson solitons were investigated . 
In that case the boundary conditions for the vector meson fluctuations had to undergo a 
special gauge transformation which corresponds to the transformation of the basis quark 
spinors described here (|4.4|). 

Before discussing the numerical treatment of eqn (|4.9|) we would like to make the 



remark that substituting the transformation (r) = ^]t(e)§^(r) into the original equa- 



tion of motion ( p.l5| ) does not result in the relation (|4.8| ) but rather yields the constraint 



= E ( ^MSgn(e^) - -erfc 



^M(^)75T-f§^(r). 



(4.15) 



Thus the equation of motion 



cannot be obtained by transforming the states but 



only by employing the Dirac equation in the rotated frame 



to extremize the energy 
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functional. The constraint ( [4 .151) does not represent an over-determination of the system 



since infinitely many states are involved. 

Due to the transcendent character of the equation of motion in the rotated frame 
(|]9D solut ions cannot be obtained for arbitrary values of the radial functions Al, At and 



V. E.g. for large distances G — requires \Al\ < \At\ in order to find a solution to eqn 
( [4.9|) . Thus the treatment of the NJL soliton in the chirally rotated frame is not very 
well suited for an iterative procedure to find the self-consistent solution. The reason is 
that a small deviation of the radial functions A^, At and V from those corresponding to 
this solution can render eqn ( [4.9| ) indissoluble for 9(r). Then it is not unexpected that at 
large distances the solution to the rotated equation of motion ( ^^91 ) becomes unstable and 
the original profile function cannot be reproduced for r > 2fm. For smaller values of r the 
original profile function is well reproduced. In figure 4.2 this behavior is displayed. The 
self-consistent solution to eqn ( p.l5| ) serves as ingredient to evaluate the radial functions 
At, At and V . The solution to eqn ( [4.9| ) is then constructed and compared to the original 
profile function. 

5. The soliton without axial vector mesons 

In this section we will examine the soliton solutions when the chirally rotated axial 
vector field is set to zero, = 0. In Skyrme type models this approach has been 
frequently used to study the structure of baryons 0, ^ |3^ . Here will ignore the effects 



of the 00 meson. There are two reasons for doing so. First, as this field represents 
an isosinglet it is not affected by the chiral rotation. The second reason is of more 
technical nature. Since the corresponding grand spin zero ansatz introduces a non- 
vanishing temporal component of a vector field the Dirac Hamiltonian in Euclidean space 
is no longer Hermitian. A stringent derivation of the corresponding Minkowski energy 
functional has unfortunately not yet been performed successfully. Several attempts for 
motivating such an energy functional have, however, been made [0, ^ ^ . Excluding 
the o^-meson we have besides the hedgehog ansatz for the chiral angle ( |1.11D only the 
Wu-Yang form for the chirally rotated vector field 



wherein 



Gir) 

and v"^ = -^emkaf-k- (5.2) 



Here G{r) refers to the dynamical p-meson field and should not be confused with the 
induced vector field discussed in eqn ( |1.13| ). From the special form of the Dirac operator 
in the chirally rotated frame ( |2.2CI| ) it is then obvious that the chiral field only appears 
in the purely mesonic part of the energy functional 

E^ = -^ Jdr |(G'(r) + 1 - cosQ{r)f + {e'{r)f + sin2e(r)| . (5.3) 



The boundary value G(r = 0) = — vr then implies that G{r = 0) = —2 p4| , |35| . The sta- 
tionary condition for the chiral angle now becomes a second order non-linear differential 
equation 

= + ;:^sin20(r) + -sine(r) (G(r) + 1 - cose(r)) . (5.4) 
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Once a profile function G{r) is provided this equation can be solved analogously to Skyrme 
model calculations. 

Undoing the chiral rotation after having set = yields a Dirac Hamiltonian which 
contains all possible grand spin zero operators of positive parity 

H = a ■ p + m/5 (cos0 + ij^r ■ rsinB) + ^^'^ ' ' ''" i^-^) 

sin© 1 

{G + 1) {t ■ (T — r ■ err ■ t) -\ {GcosQ — 1 + cosO) ck ■ (f x r). 



2r ' " ' 2r 

In account of the above mentioned boundary values for B and G this operator does 
not contain any coordinate singularities. On the contrary, the chirally rotated Dirac 
Hamiltonian becomes as simple as 

G(r) 

']-Cr = a-p + pm+ ^^ct- (r X t) (5.6) 

2r 



being, however, singular at r = 0. Most importantly eqn (|5.6| ) demonstrates that the 
eigenvalues of H (or TCr), e^, which enter the energy functional, are indeed independent 
of the chiral angle 0. Then also the contribution of the fermion determinant to the energy 
functional E'^'^^ = + i?™^ does not dependent on the chiral angle 0, i.e E'^'^^[Q,G] = 
E'^'^^[G]. Here E'^'^^ is understood as the sum of the eigenvalues in the form displayed in 
eqns (|3.5|j3.1(]| ). In particular an infinitesimal change of leaves E'^^^ unchanged. Thus 
the fermion determinant does not contribute to the equation of motion for ( |5.4| ). In 
view of the singular character of Tin this operator cannot be treated using the standard 
basis EO] but rather by employing techniques which are analogous to those developed 



in section 4. This corresponds to use the form ( ^.5|) with substituted by a reference 
profile in order to compute E'^'^^[G]. Fortunately, as the only relevant information 
concerns the boundary values of rather than its explicit form, we may equally well 
employ (|5.5|) with being the solution of ( ^.4|) . 

As the fermion determinant is a functional of G it contributes to the associated 
equation of motion 

G(r) = cos0(r) - 1 - ^| dQY.Hir)^ " x ^) ^.Mf'ie,). (5.7) 
where and ^^{r) are the eigenvalues and -functions of Hr. Furthermore 

/'(e^) = ?//.sgn(e^) - Jsgn(e^)erfc ( 



{5.t 



denotes the derivative of the fermion determinant with respect to e^. 

Before discussing the emergence of self-consistent solutions we will present a compu- 
tation of the fermion determinant associated with the Dirac Hamiltonian Tin for a give 
profile function G{r). A physically motivated profile is given by the local approximation! 
(l.a.) to the equation of motion ( [5.7[ ) (see also eqn ( |1.13| )) 

Gi-^-(r) = cos0s.c.(r) - 1 (5.9) 



°We adopt this notation from the analogous approximation in the Skyrme model which relates the 
Skyrme term to the field strength tensor of the p-meson. 
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Table 5.1. Contributions to the energy for the local approximation (5.9) as functions of the constituent 
quark mass m. All numbers are in MeV. 



m 



Cval 



^det[(^l.a.- 



rpLa.. 
Tpl.a. 



400 


500 


600 


-399 


-500 


-599 


150 


84 


61 


1221 


1513 


1480 


1372 


1597 


1541 



wherein Qs.c. i'^) refers to the self-consistent solution to the model with pseudoscalar fields 
only, which is described in section 3. In the local approximation the mesonic part of the 
energy is given by 

El:- = ^Jdr [r' {QL{r)f + 2sin^e..,.(r)} • (5.10) 

The energy eigenvalues are obtained by diagonalizing Ti with being substituted with 
the auxiliary field (f) as given in eqn ( [4 .51) . For the numbers listed in table 5.1 we have 
confirmed that these are independent of the parameter t. 

Obviously the valence quark orbit is extremely bound, its energy eigenvalue being 
approximately — m§. Numerically we observe that the upper component of the corre- 
sponding wave-function vanishes. Thus the valence quark orbit carries all features of 
an antiquark state. Nevertheless, the polarization of the Dirac sea only gives a minor 
contribution to the energy. This can be understood by noticing that a shift of the valence 
quark level from evai ~ m to Cvai ~ —ni does not change the vacuum part of the energu 
where all energy levels enter with their module, see eqn ( 3.10|) . Some caution has to be 



taken when considering the local approximation. It may be far off the actual solution 
because no local minimum of E]:^' exists. Applying Derek's theorem the chiral angle can 
be shown to collapse when G is equated to its local approximation. Nevertheless from 
the small contribution of the Dirac sea to the total energy we deduce that the fermion 
determinant only provides a small repulsive force. The local approximation furthermore 
suggests that this repulsive force decreases with increasing constituent quark mass m. 
This force should cause G to deviate from the local approximation according to the equa- 
tion of motion ( ^.71) . This deviation should in turn stablize E^ yielding a solution to 
the differential equation ( ^.4[ ). Stated otherwise: A significant deviation from the local 
approximation is needed to obtain stable self-consistent solutions. 

These features are actually reflected by the self-consistent solutions to eqns ( |5.4| , |5.7|) . 



For the specific case of m = 400MeV the profile functions are plotted in figure 5.1. Tech- 
nically we have produced this solution by diagonalizing the unrotated Dirac Hamiltonian 
( p.5| ). The resulting eigenfunctions have been transformed by f2(B) and substituted into 
the equation of motion for G (|5.7| ). We have then verified that the results are indeed 
independent of G. I.e. the calculation of the fermion determinant has been repeated 
using -|- 50 without altering the results. 50 has been chosen to satisfy the appropriate 
boundary conditions. The polarization of the Dirac sea indeed causes a sufficient devia- 
tion of G from its local approximation to yield stable solutions. On the other hand only 



6 



Here it should be noted that, as a consequence of discretizing the momentum space, the smallest 
module of an energy eigenvalue for a free quark is \J M"^ + (j:/ D)"^. The eigenvalues of free quarks are 
discussed in the appendix. 



18 



Table 5.2. Contributions to the energy for self-consistent solution as functions of the constituent quark 
mass m. All numbers are in MeV. 



m 


400 


500 


600 


Cval 


-354 


-436 


-531 




242 


170 


131 


Em 


165 


189 


157 


Etot 


407 


359 


288 



a small repulsive effect is obtained resulting in profile functions which have a very small 
spatial extension. At r ?s l/4fm the chiral angle has already dropped to 1/10 of its value 
at r = 0. In the same way the total collapse is avoided, the upper component of the 
valence quark wave-function becomes non-vanishing. Simultaneously the corresponding 
energy eigenvalue increases from —m as shown in table 5.2. I.e. particle characteristics 
are admixed. Furthermore table 5.2 shows that the part of the energy which is due to 
the fermion determinant is larger than for the local approximation (c/. table 5.1). This 
suggests that this part of the energy does not vanish for field configurations where the 
profiles are collapsing to a (5-like shape. In this way the total collapse is avoided and 
stable solutions do exist. Nevertheless the valence quark orbit remains strongly bound 
and possesses properties commonly attached to presence of the ai meson. Thus we con- 
clude that indeed the chiral invariant elimination of the ai meson carries over information 
which would be lost if this field were simply set to zero. A further feature which can 
be asserted to the presence of the ai meson is represented by the small extension of the 
soliton profiles (see figure 5.1). Previously it has been demonstrated that a non-vanishing 
A^-field provides a squeeze of the chiral angle. This can e.g. be inferred from figure 1 of 
ref.[^. Unfortunately the total energy gets very small (~ 400MeV). Hence one has to 
wonder whether this model can successfully applied to the description of baryons without 
amendments. It is obvious that a strong repulsion is called for. 

From the phenomenology of the nucleon-nucleon interaction is seems likely that the 
inclusion of the isoscalar vector meson uj may provide (at least some) repulsion. However, 
as already remarked, for the non-perturbative treatment of the NJL soliton a stringent 
derivation of the associated Minkowski space energy functional is not available. For- 
tunately the current approach allows one to incorporate an approximation to the full 
treatment of the cu-meson. Although this approximation may be somewhat crude it 
should at least provide a reliable answer to the question: Does the valence quark energy 
remain negative? 

The approximation we are going to consider relies on a power expansion of the NJL 
action in the uj- fieldi. The leading order is just the coupling to the baryon current 
Ncu^iiB^. The crudity of the approximation consists of ignoring all other terms in the 
expansion of the fermion determinant and assuming the leading order gradient expression 
for B^. The latter then becomes the topological current 

B,{U) = -i^e,.,.tr {(t/ts'^f/) (f/t^^t/) (f/t^-f/)! . (5.11) 

The mesonic part of the action. Am, contains a term quadratic in u. To this end the 

'^By construction, this is an analytical functional of to. Thus the energy is well defined in both 
Euclidean and Minkowski spaces. 
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cu-dependent parts of the Lagrangian collect up to 



li^ = 7^^,^^ + Ncuj,B^{U). (5.12) 

/(jr2 

This expression does not contain any derivative of the a;-field. Hence it may be eliminated 
by its stationary condition yielding the ordinary 6**^ order (in gradients of the chiral field) 
term 

Ce = ^elB,{U)B^{U) with el= J^'^^ , (5.13) 

Adding / dl^xC^ to the mesonic part of the action gives for the corresponding part of the 
static energy 

eg f dr f d&{r)^ ^ 



= (H) + ^ [-^ j sin^e(r). (5.14) 

In this case a stable minimum of exists even if the local approximation were assumed 
for G{r). The equation of motion for the chiral angle now becomes 



2rfe 1 . 2 . ^ 

— - = — + — sm2e + -sme (G + 1 - cosB) 

dr^ r dr r 

el rfBsin^e /sine dQ \ 

H TTTT-^ :r- — cosG . (5.15) 

27r4/^ dr \ r dr J ^ ^ 



el /sin9\^' 
An^f^ \ r J 



The solution to this equation together with (|577|) is plotted in figure 5.2. As expected the 
chiral angle receives a sizable extension. This is also reflected by a large mesonic part of 
the energy Em = 1607MeV for m = 400MeV. In the same manner the p-meson profile 
gets wider and deviates only slightly from its local approximation. As the profiles get 
extended the valence quark energy tends toward — m and regains its antiparticle character 
which was already observed in the local approximation. We conjecture that this property 
is common to all models which permit a sizable extension of G{r). Also the contribution 
of the fermion determinant to the total energy decreases. E.g. for m = 400MeV we find 
E^^* = 132MeV. Together with E^ this adds up to Etot = 1739MeV. 

Let us finally compare the results obtained in the present model (where the axial 
vector degree of freedom is eliminated in a chirally invariant way) with other treatments 
of the isovector (axial) vector mesons in the context of the NJL model. The numbers for 
the associated energies are displayed in table 5.3. 

It is obvious that either setting the original axial vector degree of freedom to zero 
{An = 0) or the chirally rotated one (A^ = 0) amounts to totally different approaches 
to the soliton sector of the NJL model. While the first one yields a positive valence 
quark energy (for a reasonable choice of parameters), the orbit is much more strongly 
bound in the second approach. This certainly represents an effect due to the axial vector 
field Afj^ ^ 0. Unfortunately the total energy of the chirally rotated treatment cannot be 
assigned to any of the other treatments due to its smallness. 

6. Conclusions 
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Table 5.3. Contributions to the energy for self-consistent solution in various treatments of the NJL 
model. Those meson fields which are allowed to be space dependent are indicated. The constituent 
quark mass m =400MeV is common. All numbers are in MeV. 



TT — p [|T^ n — p — ai [|T7l This model 





313 


-222 


-351 


^det 


711 


543 


240 


Em 


149 


393 


175 


Etot 


861 


937 


415 



We have investigated the role of chiral transformations for the evaluation of fermion 
determinants. When these transformations are topologically trivial they provide a use- 
ful tool to evaluate the chiral anomaly. Furthermore they can be used to demonstrate 
the equivalence between the hidden gauge and massive Yang-Mills approaches to vector 
mesons. As a further application of the chiral rotation we have shown that the axial 
vector degree of freedom can be eliminated without violating chiral symmetry. A gen- 
eralization to the case when chiral fields have a topological charge different from zero is 
not straightforward. Even though the special transformation we have been considering 
is unitary its topological character prevents the eigenvalues and -vectors of the original 
Dirac Hamiltonian to be regained from the rotated Hamiltonian unless the boundary 
conditions are transformed to the topologically non-trivial sector accordingly. Further- 
more we have observed that the stationary conditions to the static energy functional in 
the topologically distinct sectors are not related by the transformation of the equation 
of motion. The boundary conditions for the chiral field obtained from the stationary 
condition have been found being invariant under the chiral rotation only when the basis 
quark fields are taken from the topological sector associated with the chiral transforma- 
tion. Diagonalization of the rotated Dirac Hamiltonian in this basis can be reformulated 
into a problem where the induced vector fields ( |1.1CI| ) belong to the topologically trivial 
sector (c/. eqn ([4.6|)). In order to diagonalize the resulting operator ( |4.6| ) the standard 
basis ^ may be employed. These techniques have been applied to the case when 
the chirally rotated axial vector degree of freedom is absent. We have seen that in this 
case soliton solutions do exist and that these are totally different from those which are 
obtained in the NJL model with the original axial vector field being neglected. In par- 
ticular, the chirally rotated vector meson contains an important information inherited 
from the original axial vector meson: A negative valence quark energy. Unfortunately, 
the chirally rotated approach as it stands cannot be considered to be realistic. This is 
due to the instability of the mesonic part of the energy in the local approximation for 
the vector meson field. We have, however, made plausible that the incorporation of the 
tu-meson will provide an approach suitable to describe the physics of baryons. 

Furthermore we have seen that for field configurations wherein the rotated vector 
fields possess a reasonable spatial extension leads to a valence quark energy being of 
the amount —m. I.e. the valence quark has (almost) joined the negative Dirac sea; a 
picture which underlies Skyrme type models. Thus those amendments of the Skyrme 
model which drop the chirally rotated axial vector meson field [^4[ |3^, ^ gain strong 
support from the investigations in a microscopic quark model presented in this paper. 

Let us finally point to a possible application of the techniques developed in section 4 
and 5 for a different area of physics: The quark spectrum in the background field of an 
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instanton field configuration. This subject has gained some interest recently for the study 
of the partition function of QCD[0. Assuming temporal gauge, Aq = 0, the evaluation 
of this spectrum can be reduced to an eigenvalue problem for the Dirac Hamiltonian 



a ■ p — cx ■ A{x, 0:4) + Pnio 



(6.1) 



with the time coordinate X4 acting as a parameter. Here A^ = is a pure but 

singular gauge configuration. V may be chosen in hedgehog form relating color to coor- 
dinate space. An explicit expression is e.g. given in eqn (16.50) of ref. |Q. The similarity 
between the one-particle operators ( |6.1| ) and ( |5.6| ) is apparent^. In order to diagonalize 
( |6.1| ) the singularity carried by A has to be removed by a topologically non-trivial trans- 
formation. Since this transformation is different for various time slices a shifting of the 
quark levels as shown in figure 4.1 may occur along the path —00 < X4 < +00. 
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Appendix: Diagonlization of the Dirac Hamiltonian 

Technically the discretized eigenvalues of the Dirac Hamiltonian H ( ^.l|j3.3|J^ ) are 
obtained by restricting the space R3 to a spherical cavity of radius D and demanding 
certain boundary conditions at r = D. Eventually the continuum limit D ^ 00 has to be 
considered. In order to discuss pertinent boundary conditions it is necessary to describe 
the structure of the eigenstates of H. Due to the special form of the hedgehog ansatz the 
Dirac Hamiltonian commutes with the grand spin operator 

G = J + - = Z + - + - (A.l) 
2 2 2 ^ ^ 

where J labels the total spin and I the orbital angular momentum, t/2 and cr/2 de- 
note the isospin and spin operators, respectively. The eigenstates of H are then as well 
eigenstates of G. The latter are constructed by first coupling spin and orbital angular 
momentum to the total spin which is subsequently coupled with the isospin to the grand 
spin [0. The resulting states are denoted by \ljGM) with M being the projection of 



G. These states obey the selection rules 




(A.2) 



The Dirac Hamiltonian furthermore commutes with the parity operator. Thus the eigen- 
states of 7i with different parity and/or grand spin decouple. The coordinate space 



"Note, however, that A contains parity violating pieces. 
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representation of the eigenstates is finally given by 



[ /f •+^^)(r)|G + 1G + \GM) ) + 1, -/f '+^^)(r)|G - IG - \GM) ) ^^'''^ 

^(G,-) _ (^9f'-^'Kr)\G + IG + \GM) \ ( igf^-''\r)\G -IG- \GM) \ 

\ -ff^-''\r)\GG+\GM) ff^-^^\r)\GG~\GM) ^"^'^^ 

The second superscript labels the intrinsic parity Ilintr which enters the parity eigenvalue 
via n = (— 1)*^ X Ilintr- In the absence of the soliton {i.e. G = 0) the radial functions 
(^(G,+;i)^j,^^ j(G,+;i)^j.^^ g^^^ 0^.^ given by spherical Bessel functions. E.g. 



gf'+-'^\r) = Nk^l + m/E jaikr), fjf'+-''\r) = NkSgn{E)^l - m/ E jG+i{kr) (A.5) 

and all other radial functions vanishing represents a solution to 7i(0 = 0) with the energy 
eigenvalues E = ±\^W+Tn? and parity (— 1)*^. Nk is a normalization constant. 

Two distinct sets of boundary conditions have been considered in the literature. Orig- 
inally Kahana and Ripka proposed to discretize the momenta by enforcing those 
components of the Dirac spinors to vanish at the boundary which possess identical grand 
spin and orbital angular momentum, i.e. 

9f'^-^'\D) = 9f'^-''\D) = jf'~''\D) = ff'-^'\D) = 0. (A.6) 

This boundary condition has the advantage that for a given grand spin channel G only 
one set of basis momenta {fenc} is involved. These knc make the G^^ Bessel function 
vanish at the boundary {jaikncD) = 0). However, this boundary condition has (among 
others) the disadvantage that the matrix elements of flavor generators, like are not 
diagonal in momentum space. If the matrix elements of the flavor generators are not 
diagonal in the momenta a finite moment of inertia will result even in the absence of a 
chiral field . Stated otherwise, in this case the boundary conditions violate the flavor 



symmetry. This problem can be cured |]36[ by changing the boundary conditions for the 
states with Ilintr = — 1 

9f^^-''\D) = gf^^-^'KD) = gf'-^'\D) = gf'-^'\D) = (A.7) 

i.e. the upper components of the Dirac spinors always vanish at the boundary. The 
diagonalization of the Dirac Hamiltonian ( |3.1| ) with the condition ( |A.7| ) is technically 
less feasible since it involves three sets of basis momenta {knc-i}, {knc+i} and {knc} for 
a given grand spin channel. In table A.l we compare some properties of the two boundary 
conditions ( |A.6D and ( |A.7| ) in the case when no soliton is present. The first four quantities 



appearing in that table show up in various equations of motion when e.g. also (axial-) 
vector mesons are included In case such a quantity is non-zero the vacuum gives 



a spurious contribution to the associated equation of motion. For an iterative solution to 
the equations of motion this spurious contribution has to be subtracted |^ . It should be 



noted, however, that the relations listed in table A.l are all satisfied for both boundary 
conditions in the continuum limit D ^ oo. 

So far the discussion of the boundary conditions has only effected the point r = D. 
In the context of the local chiral rotation it is equally important to consider the wave- 
functions at r = 0. As already mentioned the solutions to the Dirac equation ( p.l3|) are 
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Table A.l. Properties of the two boundary conditions (A. 6) and ( |A.7| ) in the baryon number zero sector, 
/(r) represents an arbitrary radial function. 



Quantity 












Condition (|A.6 


) Condition (A.7|) 


^I/375^T ■ f ^^erfc ( 


A 


) sgn(e^) = 


yes 




yes 


"^la ■ (r X f) ^^erfc ( 


A 


) sgn(e/,) = 


no 




yeb 


^Ja • T^^erfc ( 


A 


) 


sgn( 




yes 




yes 


^Ja ■ rr ■ f ^^erfc ( 


A 


) 


sgn(e^) = 


yes 




yes 


tr(/5/(r)) = 












yes 




no 


tr(7-r/(r)) =0 












yes 




yes 


(yu Ti z/) = ioT kf, ^ k 










no 




yes 



given by spherical Bessel functions in the free case, = 0. Except of jo these vanish at 
the origin. In the case G 7^ we adopt the boundary conditions 0(0) = — nvr thus no 
singularity appears in the Dirac Hamiltonian ( pj.l3D at r = 0. Therefore the radial parts 
of the quark wave-functions may be expressed as linear combinations of the solutions to 
the free Dirac Hamiltonian. E.g. 



k 

ff'^'''\r) = ^ V,kmNkSgn{EkG)^l - m/EkG jG+iihcr) (A.8) 



where the eigenvectors V^fc[0] are obtained by diagonalizing the Dirac Hamiltonian in the 
free basis. It should be stressed that the use of the free basis is only applicable because 
the point singularity hidden in r ■ f has disappeared. If singularities show up for certain 
field configurations the basis for diagonalizing Ti has to be altered. This has been the 
central issue of section 4. There we have described that the eigenvalues and -vectors 
can be obtained by adopting a basis which is related by a local chiral transform to the 
one described in eqn ( [A.5| ) together with the appropriate boundary conditions (c/. eqn 
or ( A. 7 )). Here we will construct an alternative basis for the G = channel by 



explicitly transforming the boundary conditions. This has the advantage that there is no 
need to introduce the auxiliary field as in eqns ( ^.4] and |4.5| ). We therefore consider 
the application of the rotation ( |4.ij| ) at r = to the eigenstates. In the O"*" channel this 
leads to 

It is then obvious that a pertinent basis is given by 



-z sgn(E)^l - E/m ji(A;r)|OiOO) \ 

This is, of course, no longer a solution to the free Dirac equation. Moreover, a single 
component of this spinor does not even solve the free Klein Gordan equation. At r = D 
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the chiral rotation equals unity. Thus we demand the discretization condition ji{kD) = 
according to eqn (|A.6|) i.e. the momenta are taken from the set {fc„i}. With this basis 
we have succeeded in ehminating the singularities at r = and keeping track of the 
boundary conditions at r = D. We are thus enabled to numerically diagonalize T-Cr. 
When comparing with the eigenvalues of Ti we again encounter a form of the "infinite 
hotel story" : one state is missing in the negative part of the spectrum while an additional 
shows up in the positive part. The missing state turns out to be the one at the upper end 
of the negative Dirac sea, i.e. E ^ —m. Then it is important to note that in addition to 
the states with finite A;, the basis ( |A.10| ) together with the boundary condition jiikD) = 
also allows for the "constant state" with k = 0. In the continuum limit {D oo) this 
state is absent. Including, however, this state for finite D in the process of diagonalizing 
Tin finally renders the missing state. This is not surprising since application of the inverse 
chiral rotation Vl^ [r = 0) onto this "constant state" leads to an eigenstate of the free Dirac 
Hamiltonian with the eigenvalue —m. It should be remarked that for the free unrotated 
problem a "constant state" with eigenvalue —m is only allowed in the = 1~ channel. 
Although fl[Q) does commute with the grand spin operator its topological character 
mixes various grand spin channels via the boundary conditions. 

Accordingly the diagonalization of Tin in the = 0~ channel demands the basis 
states 

isgn{E)^l-E/mjoikr)mOO) 
^ + ji(A;r)|0|00) 



with the boundary condition ji{kD) = in order to be compatible with the Kahana- 
Ripka diagonalization of Ti. The additional "constant state" needed here corresponds 



to a state with eigenvalue +m of the free unrotated HamiltonianS. 

We have finally been able to diagonalize the chirally rotated Hamiltonian Tipi in the 
G = sector by very tricky means. It should also be kept in mind that there are now ad- 
ditional states at the upper (from O'^) and lower (from 0~) ends of the spectrum which do 
not possess "counterstates" in the (7 = part of the spectrum in the unrotated problem. 
For the dynamics of the problem they are of no importance because their contribution 
to physical quantities is damped by the regularization. However, their existence reflects 
the topological character of the chiral rotation. 

In the other channels i.e. G > 1 we have unfortunately not been able to construct 
a set of basis states which rendered the eigenvalues of the unrotated Hamiltonian along 
the approach described above. In the G = Q sector we have already seen that a "global 
rotation" —ir ■ r^j^ is needed for the basis states in order to accommodate the boundary 
conditions at r = 0. Furthermore a mixture of different grand spin channels appears via 
the boundary condition since this "global rotation" deviates from unity at r = D. 



"The additional "constant states" thus do not aher the trace of the Hamiltonian. 
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FIGURE CAPTIONS 



Fig. 4.1 

A schematic plot of the spectrum of the rotated Hamiltonian T-CnR = fi(n0)7Yfi^(?20). 
Fig. 4.2 

Comparison of the self-consistent profile in the unrotated formulation (dashed line) and 
the solution to eqn (|4.9|) (solid line). 
Fig. 5.1 

The meson profiles (left) and valence quark wave-function (right) of the self-consistent 
solution to eqns (|5.4| , ^.7|) . Also shown is the local approximation to the vector meson 
field G^'^' = cosG — 1. Here the constituent quark mass m =400MeV is assumed. 
Fig. 5.2 

Same as figure 5.1 with the 6*'^-order term ( ^.14|) included. The upper component of the 
valence quark wave-function is negligibly small. 
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